ASYMPTOTICALLY EFFICIENT TRIANGULATIONS OF THE <i-CUBE 



DAVID ORDEN AND FRANCISCO SANTOS 



Abstract. Let P and Q be polytopes, the first of "low" dimension and the second of 
"high" dimension. We show how to triangulate the product PxQ efficiently (i.e., with few 
simplices) starting with a given triangulation of Q. Our method has a computational 
part, where we need to compute an efficient triangulation of P X A m , for a (small) 
natural number m of our choice. A m denotes the m-simplex. Our procedure can be 
applied to obtain (asymptotically) efficient triangulations of the cube I": We decompose 
jn _ jk x jn—k^ £ or a sma n Then we recursively assume we have obtained an efficient 
triangulation of the second factor and use our method to triangulate the product. The 
outcome is that using k = 3 and m = 2, we can triangulate I n with O(0.816 n n!) simplices, 
instead of the O(0.840 n n!) achievable before. 

Keywords: dissection, triangulation, size, cube, efficiency, simplcxity. 



1. Introduction 

"Simple" triangulations of the regular c£-cube I d — [0, l] d have several applications, such 
as solving differential equations by finite element methods or calculating fixed points. See, 
for example, In particular, it has brought special attention both from a theoretical 

point of view and from an applied one to determine the smallest size of a triangulation of 
the d-cube (see Section 14.5.2] for a recent survey). Let us point up that the general 
problem of computing the smallest triangulation of an arbitrary polytope is NP-complete 
even when restricted to dimension 3, see j2]. 

When we speak about triangulations of a polytope P of dimension al we mean decompo- 
sitions of P into ci-simplices that (i) use as vertices only vertices of P, and (ii) intersect face 
to face (i.e., forming a geometric simplicial complex). Some authors do not require these 
two conditions in triangulations. We will always require the first one, and when the second 
condition is not fulfilled, we call the decompositions simplicial dissections of P. The size 
of a triangulation or dissection T is its number of d-simplices and we denote it \T\. It is 
an open question whether high dimensional cubes admit dissections with less simplices than 
needed in a triangulation. Actually, the minimum size of dissections of I is unknown, while 
the minimum triangulation is known (see below). 

The paper |2j describes a general method to obtain the smallest triangulation of a polytope 
P as the optimal integer solution of a certain linear program. The linear program has as 
many variables as d-simplices with vertex set contained in the vertices of P exist. That is, 
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( dimfpyt-f) ^ ^ e ver tices of P are in general position and less than that if not. For the 
d-cube, the direct application of this method is impossible in practice beyond dimension 
4 or 5. With a somewhat similar method but simplifying the system of equations via the 
symmetries of the cube, Anderson and Hughes 1 have calculated the smallest size among 
triangulations of the 6-cube and the 7-cube, in a computational tour-de-force which involved 
a problem with 1 456 318 variables and ad hoc ways of decomposing the system into smaller 
subsystems. The smallest sizes up to dimension 7 is shown in Table 

In order to compare sizes of triangulations of cubes in different dimensions, Todd [TO] 
defines the efficiency of a triangulation T of the d-cube to be the number (\T\/d\)~x. This 
number is at most 1, since every simplex with integer vertices has a multiple of l/d\ as 
Euclidean volume, so |T| < dl. Triangulations of efficiency 1 (i.e., unimodular) can be 
easily constructed in any dimension. On the other extreme, Hadamard's inequality for 
determinants of matrices with coefficients in [—1, 1] implies that the volume of every d- 
simplex inscribed in the regular d-cube I d = [0, l] d is at most (d + l)^ d+1 ^ 2 /2 d d\. Hence, 
every triangulation has size at least 2 d d\ / (d+ lY d+1 ^ 2 and efficiency at least 2/ '(d+l)~32~ « 

2/v^TT. 

Following the notation in [HI Section 14.5.2], let <f>d and pd be the smallest size and 
efficiency, respectively, of all triangulations of the cube of dimension d. The number <f>d (or 
some variations in which one or both of the conditions (i) and (ii) are not required) is known 
as simplexity of the d-cube. Obviously, pd — (0 c ;/d!) 1 / d . 

In [H] (see also p], pages 283-284), Haiman observes that a triangulation of I k+l with 
tkti ( k ^ 1 ) simplices can be constructed from given triangulations of I k and with tk and t\ 
simplices respectively. With this, one easily concludes: 

Theorem 1.1 (Haiman). For every k and I, pk+i k+l < Pk k pi ■ 

Corollary 1.2. The sequence (pi) e pj[ converges and 

Urn pi < p d VdE N. 

Proof. Let us fix d € N and k G {1, . . . , d}. Haiman's theorem implies that, for every i E N, 

p k+ld <p k k/{k+id) Pd u/{k+id) . 

Since the right-hand side converges to p d when i grows, the d subsequences of indices modulo 
d, and hence the whole sequence (pi) . pj, have upper limit bounded by pd- An upper limit 
bounded by every term in the sequence must coincide with the lower limit. □ 

What is the limit of this sequence? In particular, is it positive or is it zero? The known 
values of pd (up to d = 7) form a strictly decreasing sequence, as shown in Table ^ but it 
is not even known whether this occurs in general. 

Concerning lower bounds, the only significant improvement to Hadamard's inequality 
has been obtained in [S], where the same volume argument is used, but with respect to a 
hyperbolic metric. The last row of the Tabled shows the lower bound obtained, translated 
into efficiency of triangulations. For small dimensions, it is an excellent approximation of 
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Dimension 


1 


2 


3 


4 


5 


6 


7 


8 


Smallest size {<j>d) 


1 


2 


5 


16 


67 


308 


1493 


< 11944 


Smallest efficiency (p,i) 


1 


1 


.941 


.904 


.890 


.868 


.840 


< .859 


Lower bound Hadamard's ineq. 


1 


.877 


.794 


.731 


.683 


.643 


.610 


.581 


Lower bound Smith 


1 


1 


.941 


.889 


.833 


.789 


.751 


.718 



Table 1. Smallest size and efficiency of triangulations of I d . 



the smallest efficiency. Asymptotically, it only increases the lower bound obtained using 
Hadamard's inequality by a constant factor y/Z/2. 

In this paper we propose a method to obtain efficient triangulations of a product polytope 
P x Q starting from a triangulation of Q and another of P x A m_1 , where A m_1 denotes 
a simplex of dimension m — 1 and m is any relatively small number. We apply this with 
P being a small-dimensional cube and Q a high-dimensional one, iteratively. This allows 
to obtain asymptotically efficient triangulations of arbitrarily high-dimensional cubes from 
any (efficient) triangulation of I 1 x A m . Sections 3, 4, and 5 explain our method, which 
is first outlined in Section 2. If the reader is happy with dissections, we cannot offer better 
efficiencies for them than for triangulations, but at least he or she can skip Section 5, which 
contains most of the technicalities in this paper. 

The asymptotic efficiency of the triangulations obtained, clearly, depends on how good 
our triangulation of P x A m_1 is. Finding the triangulation of P x A" 1-1 which is optimal 
for our purposes reduces to an integer programming problem, similar to finding the smallest 
triangulation of that polytope (actually, it is the same system of linear equations, with a 
different objective function). Using the linear programming software CPLEX, we have solved 
the system for some values of I and m. The best triangulation we have found is one of 
I x A 2 , with which we obtain 

lim Pl < \ » 0.8159 

The best bound existing before was limj_,oo pi < pi = 0.840. In other words, we prove 
that (asymptotically) the d-cube can be triangulated with 0.8159 d d! simplices, instead of 
the 0.840 d d! achievable before. 

It has to be observed that, even if the particular triangulation of I 3 x A 2 that we use 
was obtained by an intensive computer calculation, once the triangulation is found it is a 
simple task to check that it is indeed a triangulation and compute the asymptotic efficiency 
obtained from it. Actually, in Section [7| we use the so-called Cay ley Trick [7] to do this 
checking with no need of computers at all. We also use the Cayley Trick to explore the 
minimum efficiency that can be obtained from the product I 2 x A fe for any k. We briefly 
explain the trick in Sectional where we also interpret our whole construction in terms of it. 

2. Overview of the method and results 

We start describing Haiman's proof of Theorem ll.il which is related to our method. Let 
Tfe and TJ be triangulations of the regular cubes I k and P , respectively. The product Tk x 1} 
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of the two triangulations gives a decomposition of the cube I k x F = I k+l into |Tfc| • |T;| 
subpolytopes, each of them isomorphic to the product of simplices A k x A 1 . It is well known 
that every triangulation of A fc x A' has size ( k ^ 1 ) (see, e.g., Chapter 7]). Hence, refining 
Tfc x Ti in an arbitrary way one gets a triangulation of I k+l of size |Tfc| • \T{\ ■ ( t ). The 
implication of this is that starting with a triangulation of a cube I k of a certain efficiency p, 
one can construct a sequence of triangulations of I nk for n € M whose asymptotic efficiencies 
converge to p. 

Our method is, in a way, similar. Starting from a triangulation of J™ -1 and another one 
of x A m_1 , we get one of with the following general method to triangulate P x Q 

starting from a triangulation of Q and another one of P x A m_1 : 

(1) We first show (Sections and 0} how to obtain triangulations of P x A" -1 from 
triangulations of P x A" 1-1 , where n — 1 = dim(Q) is supposed to be much bigger 
than m — 1. We call our triangulations multi-staircase triangulations. 

(2) A triangulation of Q induces, as in Haiman's method, a decomposition of P x Q 
into polytopes isomorphic to P x A™ . Each of them can be triangulated using 
the previous paragraph, although this in principle only gives a dissection of P x Q; 
if we want a triangulation, we have to apply (1) to all the polytopes P x A™ -1 in 
a compatible way. In Section [S] we will show how to do this using an m-coloring of 
the vertices of Q. 

(3) The analysis of the size of the triangulation obtained will also be carried out in 
Section CD 

In step (2), the final size of the triangulation is just the sum of the individual sizes of the 
triangulations used for the different subpolytopes P x A" -1 . In particular, if we are just 
interested in obtaining dissections, we can take an efficient triangulation of P x A™ -1 and 
repeat it in every subpolytope. 

In step (1) the computation of the size is more complicated and to state it in a simple 
way we introduce the following definitions. When we speak about simplices in a polytope P 
we will implicitly suppose that its vertices are vertices of P, and we will identify the simplex 
with its vertex set. 

Definition 2.1. Let P be a polytope of dimension I and let r be a simplex of dimension 
I + m — 1 in P X A" 1-1 . Let {vi, . . . , v m } be the vertices of A m_1 . Then, r, understood as 
a vertex set, decomposes as r = t\ U • • • U r m with TjCPx {vi}- 

(i) We define the sequence (|ti| — 1, ... , \r m \ — 1) to be the type of the simplex r. The 
weight of a simplex r of type (ti, . . . , t m ) is the number jy^— ■ 

(ii) The weighted size of a triangulation T of P x A" 1 ^ 1 is J^reT weight (r), and the 
weighted efficiency is 

Ergr wei g ht ( T ) 
m l 

With this, our main result can be stated as: 

Theorem 2.2. Consider polytopes P of dimension I and Q of dimension n — 1. Let m 
be such that m < n. Given a triangulation Tq of P x A m_1 of weighted size to and a 
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triangulation Tq of Q, then there are triangulations of P x Q with size at most 

Using this method to triangulate p+ n ^ 1 = p x J™ -1 starting from a triangulation of 
I" -1 and another one of P x A m_1 we will conclude: 

Theorem 2.3. If there is a triangulation To of P x A m_1 with weighted efficiency po, then 
for every e > and all n bigger than Impo/e we have 

Pn+l-i n+l - X <Pn-i n -\po + e) 1 . 

As a consequence, 

lim pt< po- 

i — >00 

Observe that the definition of weighted efficiency makes sense even in the case m = 1, 
where all the simplices of P x A have the same type, equal to (I) if dim(P) = I, and the 
same weight, equal to l/l\. In particular, the weighted efficiency of a triangulation of x A 
is the same as the usual "non-weighted" one. Another common point between efficiency and 
weighted efficiency is that the weighted efficiency of a triangulation of / x A" 1-1 is always 
less or equal to 1, and it is 1 if and only if the triangulation is unimodular; i.e., if every 
simplex has volume 1/(1 + m — 1)! (this is proved in Section EJ). 

We now describe the practical results obtained. The last theorem leads us to study the 
smallest weighted efficiency of triangulations of x A" 1-1 ; let us denote it by p;, m . This 
number can be calculated minimizing a linear form over the so-called universal polytope of 
all the triangulations of x A m_1 . 

The definition of this universal polytope for triangulations of an arbitrary polytope P is 
as follows (see jH) for further details): let S(P) be the set of all the simplices of maximal 
dimension which use as vertices only vertices of P. Given a triangulation T of P, its incidence 
vector Vt & has a 1 in the coordinates corresponding to simplices of T and a in the 

others. The universal polytope of P is conv{Vr : T is a triangulation of P}. 

In our case, to calculate the smallest weighted size (and efficiency), we have to minimize 
over the universal polytope of P x A m_1 the linear form having as coefficient of each sim- 
plex its weight l/rii=i (To minimize non-weighted efficiency and size one just uses the 
functional with all coefficients equal to 1). One of the results in [H] is a description of the 
vertices of the universal polytope as integer solutions of a certain system of linear inequali- 
ties derived from the oriented matroid of P. Thus, our minimization problem is restated as 
an integer linear programming problem. 

In order to apply our method, we have used the program UNIVERSAL BUILDER by Jesus A. 
de Loera and Samuel Peterson 0]. Given as input the vertices of P, the program generates 
the linear system of equations defining the universal polytope of P. The output is a file 
readable by the linear programming software CPLEX. We have created a routine that generates 
our particular objective function. Table [5] shows the results obtained in the cases we have 
been able to solve. Note that pi t i = pi and so the column m = 1 in Table is the same as 
the second row of Table H 
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The fact that p\, m — 1 for every m reflects that every triangulation of the prism / x A m 1 
is unimodular. In the case of P2,m we prove (Subsection 17.1(1 that the smallest weighted 
efficiency is always a/|~3t7i 2 /4] /m 2 . That is to say, -\/3/4 for even m and -\/3/4 + 0(m~ 2 ) 
for odd m. 

The computation of ^3,3 involved a system with 74400 variables, whose resolution by 
CPLEX took 37 hours of CPU on a SUN UltraSparc. 



I \ m 


1 


2 


3 


> 3 


1 


1 


1 


1 


1 


2 


1 


y| re 0.866 


^| w 0.8819 




3 


{J^ w 0.941 


^ 8355 


^^Z» a 0.8159 





Table 2. Values of /5; iTO for / < 2 or for I = 3 and m < 3. 



3. Polyhedral subdivision of P x A fcl+ "'+' Cm 1 induced by a polyhedral 

SUBDIVISION OF P X A m_1 

We call polyhedral subdivisions of a polytope P its face-to-face partitions into subpolytopes 
which only use vertices of P as vertices. 

Let P be a polytope of dimension I. Let m and ki, . . . , fc m be natural numbers and let 
us call n := k\ + ■ ■ ■ + k m . Let Uj., . . . ,v m be the m vertices of the standard simplex A™ -1 

and v{, . . . , v\ , , v™ , ■ ■ ■ , v™ the vertices of A" -1 . Observe that, implicitly, we have 

the following surjective map: 

vert(A™- 1 ) -> vert(A m - 1 ) 

Vj H-» I7j 

This map uniquely extends to an affine projection ttq : A' 1-1 — > A m , In turn, this induces 
a projection 

7T = 1 x 7T : P x A"" 1 -> P x A™" 1 
(p,o) i-> (p,7r (a)) 

Given the projection 7r and a polyhedral subdivision S of the target polytope P x A m_1 , 
it is obvious that the inverse images 7r _1 (i?) of the subpolytopes of S form a subdivision of 
P x A™" 1 into subpolytopes matching face to face. In a more general projection it would 
not be true that those subpolytopes use as vertices only vertices of P x A n . But it is true 
in our case: 

Lemma 3.1. Let B C P x A" 1-1 be a subpolytope with vert(P) C vert(P x A m_1 ). Let n 
be the projection considered before. Let be B — {(p, Vj) : (p,Vi) € vert(B)}. Then: 

7r _1 (B) C conv(i3). 

Proof. Let (p,a) be a point of 7r _1 (P), so (p, 7To(a)) G B. Let us write a as a convex 
combination of the vertices of A™" 1 , that is a = Y^JiLi 2j=i ^ l j v j' with A* > 0, Vi, j and 



ASYMPTOTICALLY EFFICIENT TRIANGULATIONS OF THE d-CVBE 



7 



J2i J2j = !• Without loss of generality, we can suppose that none of the sums YljLi ^ s 
zero, otherwise everything "happens" on a face of A m_1 and we can restrict the statement 
to that face. 

Let (p\,vi), . . . , (pj^vi), (p\,v 2 ), . ■ • , (pf 2 ,v 2 ), ■ ■ ., (pT ,v m ), ■ ■ • , ipT m ,v m ) be the vertices 
of B, so we have B = {{p l hl v'j) : i — l,...,m,j = 1, . . . , k i: h = 1, . . . , li}. We write now 
(p, 7To(a)) as convex combination of the vertices of B, that is, 

m li 

(p,7T (a)) = J2J2^h (Ph> V i)' 
i=l h=l 

with fi\ > 0, Vi,h, and YliY^hV-h = 1 - Observe that Y,jLi Xj = S/i=i A 1 ?! for every i, 
because 7r (a) = Y^iLx Y^j=i ^j v i- Then, it is easy to check that: 

m ki h \i i m kt U -yi i 

i=l j=l h=l 2-rj=l A J i=l j=l ft=l Zjj=1 A j 

That is, (p, a) is a convex combination of points of B, as we wanted to prove. □ 
Corollary 3.2. Under the previous conditions; 

(i) n~ 1 (B) = conv(7r _1 (vcrt(_B))). 

(ii) vert(7r- x (S)) = B. 

Proof. In both equalities, the inclusion from right to left is trivial. In the second one, observe 
that if (p, Vj) is in B, then (p,Vi) is a vertex of B. Therefore, (p,Vj) 6 tt~ 1 (B) and it is a 
vertex of P x A™ -1 , which implies that it is a vertex of 7r _1 (.B). 
Inclusions from left to right follow from Lemma II because: 

(1) B C Tr-^vert^)) =*> Tr" 1 ^) C conv(P) C conv(7r- 1 (vert(B))). 

(2) n~ 1 (B) C conv(P) vert(7r _1 (P)) C vert(conv(P)) = B, where the last equality 
follows from the fact that all the elements of B are vertices of P x A" -1 . □ 

Corollary 3.3. Every polyhedral subdivision S of Px A m_1 induces a polyhedral subdivision 
S := 7r _1 (5) = {ir~ l {B) : B 6 S} o/Px A™" 1 . Furthermore, the vertices of each subpolytope 
7T -1 (.B) in this subdivision are B :— {(p, uj) : (p,Vi) € vert(P)}. □ 

4. TRIANGULATION OF P X A fclH ^ fc ™-! INDUCED BY A TRIANGULATION OF P X A™ 1-1 

We will suppose now that the polyhedral subdivision S of P x A™ 1-1 is a triangulation. 
A convenient graphic representation of the vertices of the polytope P x A m ~ 1 is as a grid 
whose rows represent vertices of P and whose m columns represent the vertices vi , . . . , v m of 
A m_1 . In order to represent a subset of vertices of Px A m_1 we just mark the corresponding 
squares. In the same way we can represent the vertices of P x A felH \-K m -i^ b u t now ft j s 
convenient to divide the grid horizontally in blocks, each of them corresponding to each ki 
and containing the vertices v\,...,v\. of A felH ^ fe ™-i^ 

Figure ^ shows how to obtain, with the notation of the previous section, the set B 
associated to a simplex B G S in this graphic representation. In the i-th block, rows 
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corresponding to vertices (p, in B have all its squares marked. Restricting B to that 
block gives precisely conv({p\, . . . x conv({v\, . . . , v l k .}), with the notation used for 

the vertices of B in the proof of Lemma Hill 



B B 



FIGURE 1. How to obtain B from B. 

Since B is a simplex, any subset of its vertices forms also a simplex; thus, the restriction 
of B to each block is a product of two simplices. Let us recall that the staircase triangulation 
of the product of two simplices A fc x A 1 is the one whose ( k ^ 1 ) simplices are all the possible 
monotone staircases in a grid of size (k + 1) x (I + 1). Following this analogy, we define 
multi-staircases as follows: 

Definition 4.1. Let B = n~ 1 (B) be a subpolytope of Fx A felH \-k m -i ) of the kind obtained 
in Section |31 

(i) A multi-staircase in B is any subset of vertices which restricted to every block forms 
a monotone staircase. 

(ii) The multi- staircase triangulation of B is the one which has as simplices the different 
multi- staircases (see Figure EJ). 



Figure 2. The four multi-staircases (right) forming the multi-staircase tri- 
angulation of the polytope B in the left. 



Lemma 4.2. The multi- staircases indeed form a triangulation of B and taking the multi- 
staircase triangulations of the different B's obtained from a triangulation of P x A" 1-1 we 
get a triangulation of P x A kl ^ Msm— w hich we call multi-staircase triangulation. 
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Proof. It is clear that multi-staircases form full-dimensional simpliccs in B. A way to prove 
that a collection T of full-dimensional simplices in a polytope B is a triangulation is to show 
that: 

(1) They induce a triangulation on one face of B. 

(2) For each full-dimensional simplex in the collection, the removal of any single vertex 
produces a codimension one simplex either 

• lying on a facet of B and not contained in any other simplex of T, or 

• contained in exactly another full-dimensional simplex of T which is separated 
from the first one by their common facet. 

In our case, the first condition follows by induction on k\ +• • • + k m , the base case being 
fci H h k m = m. 

For the second condition, let a be a multi-staircase, and let (p, v) be a vertex in it. Then, 
one of the following three things happens (Figure |3| gives an example of a multi-staircase): 

• If (p, v) is the only point of a in its column, then no other multi-staircase contains 
a \ {(p, v)}. Removing (p, v) produces indeed a codimension one simplex contained 
in a facet P x £k 1+ -+k m -2 of p x ^k 1+ -+k m -i _ 

• If (p, v) is the only point of a in a row within a block, then no other multi-staircase 
contains a \ {(p, v)}. Removing (p, v) produces indeed a codimension one simplex 
in a facet of B of the form 7r -1 (-B \ {(p, tto(v ))}) (remember that B = tt(B) is a 
simplex) . 

• If (p, v) is an elbow in the multi-staircase, then removing it leads to a unique different 
way of completing the multi-staircase. More precisely, let (p' , v) and (p, v 1 ) be the 
points of the multi-staircase adjacent to (p,v). Then removing (p, v) and inserting 
(p',v') produces the other possible multi-staircase. The obvious affme dependency 
(p, v) + (p', v') = (p, v') + (p' , v) implies that the two multi-staircases lie in opposite 
sides of their common facet. 

□ 



Figure 3. An example of a multi-staircase. 

Lemma 4.3. Let U be the number of vertices of B in the i-th column. Then the multi- 
staircase triangulation of B has exactly Ol=i ( ) simplices. 

Proof. In each sub-block there are { ~u^_^~ 1 ) possible monotone staircases. □ 
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5. A TRIANGULATION OF P X Q 



We will consider polytopes P of dimension I and Q of dimension n — 1. We assume we are 
given a triangulation Tq of Q, which induces a decomposition of P x Q into cells isomorphic 
to P x A" -1 . Then, any decomposition n := fci + • ■ ■ + fc m allows us to apply the procedure 
of Sections El and 0] to triangulate the cells P x A™ -1 starting from a triangulation Tq of 



There are two important tasks remaining: First, show that the triangulations of the 
different P x A™ -1 can be achieved in a coordinated way to obtain a real triangulation of 
P x Q; second, analyze the efficiency of that triangulation. Both of them will be done in 
this section, using the following trick: 

Consider a partition of the vertices of Q into m "colors" . Then, in each subpolytope 
P x A™ -1 of P x Q the vertices of the factor A™ -1 are colored as well. We use this coloring 
to construct the projections A™ -1 — > A m_1 we need for each of them. Then, on each 
common face (isomorphic to P x A k , k < n — 1) of two of the cells P x A" -1 we get the 
multi-staircase triangulation induced by the m-coloring of Q restricted to that face. Hence, 
the triangulations of the different cells P x A™ -1 intersect face-to- face, and we obtain a 
triangulation Tp X Q of P x Q, which we call multi- staircase triangulation of P x Q. 

In order to analyze the size of Tp X Q, we will suppose that the m-coloring of the vertices 
of Q is chosen at random with a uniform distribution. 

For each a € Tq let o~i := {vertices of a colored i}. And for each r 6 To, r,; := 
{vertices of r over the i-th vertex of A™" 1 }. By Lemma 14.31 the triangulation Tp X Q we 
obtain has 



simplices. The expected value of this sum, when the coloring is random, equals the sum of 
the expected values. So let us fix a pair of simplices t £ To and a 6 Tq and calculate the 
expected value of 



We call li := |t»| and ki :— \<Ji\. The l^s are considered constants, while the ki's are random 
variables depending on the coloring. They follow a multinomial distribution, in which the 
probability of the m-tuple (k\, . . . , k m ) is P{k\, . . . , k m ) = fci ,...^, m „ . 



P x A 



m — 1 





Since (; 
can write: 



x 



) 



x—/nl, where x— : 



x{x 



1) 



■ {x — n + 1) is the n-th falling power of we 




hi 



i+h-i 

h - 1 



rr=i(^-i)< 



For the numerator in the RHS of this equation we will use the following result from [TT] 
(Theorem 4.4.4): 
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Theorem 5.1. Let X\, . . . , X m be scalar random variables. Then, the expected value of 
n«=i X- 1 - is given by the formal calculus 



in which extra variables z%, i = 1, . . . , m appear with formal purposes. 

Lemma 5.2. Let l\,...,l m be positive integers with = I + m and let k\, . 
random variables obeying a multinomial distribution with — n. Then, 



be 



£(^-1 + ^-1)^) < 



Proof. If some U equals 1 we can neglect it in the statement: we remove it and will still have 
J2(h — 1) = i. Hence we will assume U > 2 for every i. 

We will use the following formula for expectations under a multinomial distribution, again 
taken from (p. 53, last formula in the proof of Theorem 4.2.1): 



m / m \ 1 

*(IK'>= E^. 

i=l \i=l / 



Since ■ 1 (lj — 1) = I and only the ki depend on the random process, applying Theorem 
15.11 to the random variables := ki + h — 2 and using the equality above we get: 



E(Y[(h - I + h - l) 1 -^- ) = 



-E(UT=i 



^(n™i^- 2 ^(n™i^) 



(ic (e™, ±«.)") ] M < n («■-»+ s) s 



< 



no 



n 
rn 



1 + 



where the last inequality comes from ij < I + 1 (since XX'i — 1) = J and l.- L — 1 > 1), and the 
one marked with an asterisk needs to be proved. For this we call % m (z) :— YiiLi z i i_2 
and G n (z) := (Y^Li m z 0™- Using that 



3=0 



j J dxi dx k i ' 



we come up to 



" Z i 3=0 



(Z< - Vp-±± ( - J ni F,,,...^ i B (?)G n -j(J). 
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Therefore, since 



d' 



-h(z) 



d 



and 
we get 



( a' 1 - 1 



J z=l 



h-l im-l / m 



j 1= j m =0 \i=l v 7 v 7 



< 



l\— 1 ^TTl — 1 / 771 // 1 

E-E nr: ) 

ii=0 j m =0 Vi=l 
ii-1 /m-1 / m 

e -e n 

31=0 j m =0 \i=l 
'/i-I 



/ i-y. f±V )„>:, . 



ii - 1 



(i 4 - 2) 



m . 



H(E( : 'V'< ^' ' ' »" 



i=l \ji=0 



n( 

i=l 



it— l 



as wanted. 



□ 



This lemma is crucial to prove the two theorems announced in Section 2. Indeed, Theorem 
12.21 is just a version of the following more precise statement: 

Lemma 5.3. Consider polytopes P of dimension I and Q of dimension n — 1. Let m be such 
that m < n. Given a triangulation Tq of P x A m_1 of weighted size to and a triangulation 
Tq °f Q> the expected size of the multi-staircase triangulation Tp X Q of P x Q is bounded 
above by 



\T Q \to{~ + l 



Proof. Lemma 15.21 implies that 
% -l + k-l* " 



* n 



k-i 



(k-iy. 



n£iOi-i) 



V m/ 



Hence: 



v <tGTq tGT i=l 



atift-i) 



™(£nd^)(^0'-™.(£ + 



1 v m / 



□ 
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Proof of Theorem \2.3l Let to — p^m 1 be the weighted size of the triangulation in the state- 
ment. For any n > ^p-: 

m ( n \ 

Pa + e > Po— [ h I ) , 

n \m / 

n l {po + ef > pW (—+1) =to(—+l) , 



(po + e)' > to (— + I 

\771 



(n — 1)! Vto 
Theorem E21 with Q = P 1-1 and P = P, tells us that 

/ n \ l 

<Pn+l-l < Pn-lto ( hi) ■ 



<t (- + l) 



or, in other words, 

(n + l-l)\p n+l - 1 n+l - 1 

(n - l)!^-!"" 1 ~ u Vm 
So, we have proved the first part of the theorem: 

Ve > , Vn>^, p„ +i _ 1 "+ i - 1 <p„_ 1 »- 1 (p + e )'. 
e 

The second part of the statement follows from the first one with arguments similar to 
those of Corollary 1 1.21 Recursively, we have that: 

Ve>0, VzeN, Vn>^p, p n+il < p n ^fn ( po + e )s¥n . 
This implies that for the given I and any fixed n, taking e = SMo. we g e t 

hm p n+a < po H . 

That is, the sequence of indices congruent to n modulo Z has limit bounded by the right 
hand side. Since we can make n as big as we want and the rest of the right-hand side are 
constants, the I subsequences of indices modulo / have limit bounded by po, hence the limit 
of the whole sequence has this bound. □ 

6. Interpretation of our method via the Cayley Trick 

The Cayley Trick allows to study triangulations of a product P x A™ -1 as mixed subdi- 
visions of the Minkowski sum P + • • • + P (n summands). We overview here this method, 
but the reader should look at [7] for more details. 

Let Qi, . . . , Q m C M d be convex polytopes of vertex sets At- Consider their Minkowski 
sum, defined as 

TCI 



1=1 



Qi { 3*1 ~f~ " " " ~t~ ■ £ ■ 



We understand Y^iLi Qi as a marked polytope, whose associated point configuration is 
Sl=i := {Si + ' ' ' + Ira ■ Qi S Ai\. Here a marked polytope is a pair (P, A.) where P is 
a polytope and A is a finite set of points of P including all the vertices. Subdivisions of a 
marked polytope are defined in Chapter 7] (sometimes they are called subdivisions of A) . 
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Roughly speaking, they are the polyhedral subdivisions of P which use only elements of A 
as vertices (but perhaps not all of them) . They form a poset under the refinement relation. 
The minimal elements are the triangulations of A. 

A subset B of YliLi is called mixed if B = B\ -| — • + B m for some non-empty subsets 
Bi C At, i = 1, . . . ,m. A mixed subdivision of X)"=i Qi ^ s a subdivision of it whose cells 
are all mixed. Mixed subdivisions form a subposet of the poset of all subdivisions, whose 
minimal elements are called fine mixed, in which every mixed cell is fine, i.e., does not 
properly contain any other mixed cell. 

We call Cayley embedding of {Qi, . . . , Q m } the marked polytope 

(C(Q 1: Q m ),C(A 1: . . .,A m j) in R rf x R'"" 1 

defined as follows: let ei, . . . , e m denote an affine basis in R m and Hi : W -> R x R 
be the inclusion given by fJ-i(x) = {x, ei). Then we define 

C{Ai,...,A m ) := U™ C(Qi,...,Q m ) := com(C(Ai, . . . ,A m )) 

Each Qi is naturally embedded as a face in C(Q\, . . . , Q m )- Moreover, the vertex set of 
C(Qi, . . . , Q m ) is the disjoint union of the vertices of all the QiS. This induces the following 
bijection between cells in C(Qi, • ■ ■ , Q m ) and mixed cells in Q± + ■ • • + Q m - To each mixed 
cell B\ + - ■ - + B m we associate the disjoint union BiU- ■ -U£> m . To a cell B in C(Qi, . . . , Q m ), 
we associate the Minkowski sum B\ -| + B m , where Bi = B D Qi. 

Theorem 6.1 (The Cayley Trick,0). Let Qi,...,Q m C R be convex polytopes. The 
bijection just exhibited induces an isomorphism between the poset of all subdivisions of 
C(Qi, . . . ,Qm) and the poset of mixed subdivisions of X^=i Qi- ^ n isomorphism tri- 
angulations correspond to fine mixed subdivisions. 

Remark 6.2. With the previous definitions, for any polytope P: 

P x A'"" 1 =C(P,.T\,P). 

In particular, in our context the Cayley Trick provides the following bijections between tri- 
angulations and mixed subdivisions: 



Triangulation of P x A r 



Triangulation of P x A r ' 



J Cayley Trick 



\ Cayley Trick 



Fine mixed subdivision of P+ ■ ■ ■ +P 



Fine mixed subdivision of P-\ 



+P 



Our interest in the Cayley Trick is two-fold. On the one hand, it provides a way to 
visualize our candidate triangulations of / x A m_1 as objects in dimension I, instead of 
I + m — 1. We will use this in Section [7| 

On the other hand, the construction of the previous sections has a simple geometric 
interpretation in terms of the Cayley Trick. More precisely, the polyhedral subdivision of 
Corollary 13.31 can be obtained as follows: Let 5 be a polyhedral subdivision of P x A" 1-1 , 
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and Sm the corresponding mixed subdivision of P+ ■ ■ ■ +P. Each cell in S decomposes 
uniquely as 

U™ 1 B i x {«(}, 

where the Bi are subsets of vertices of P. The corresponding cell in Sm is just B\ + ■ ■ ■ + B m . 
To construct our polyhedral subdivision of P x A™ -1 we just need to scale each summand 
Bi by the integer ki which tells us how many vertices of A™ -1 correspond to the vertex Vi 
of A" 1-1 . 

That is to say, from Sm we construct a mixed subdivision of P+ ■ ™ ■ +P by the formula 

' fern 

Sat = {-E>i + • • • +B\ + + B m + ■ ■ ■ +B m : B\ -\ + B m e Sm}- 

The polyhedral subdivision S of P x A" -1 stated in Corollary 13. 31 is the one corresponding 
via the Cayley Trick to the mixed subdivision Sm of P+ ■ ■ ■ +P. 

Also the type and weight of a simplex in P x A m_1 have a simple interpretation via 
the Cayley Trick. With the notation of Definition 12.11 let r = Ty U • • • U r m be a simplex. 
The corresponding cell tm in P+ ■ ■ ■ +P is the Minkowski sum of the simplices t\ , . . . , T rn , 
which lie in complementary affine subspaces. Hence, tm is combinatorially a product of m 
simplices, of dimensions t\, . . . ,t m where (tx, . . . ,t m ) is the type of t. Then the weight of 
r represents the volume of tm, normalized with respect to the unit parallelepiped in the 
lattice spanned by the vertices of tm- With this we can prove: 

Proposition 6.3. Let P be a lattice polytope of dimension I. Let V be its volume, normalized 
to the unit parallelepiped in the lattice. Then, the weighted size of a triangulation of Px A m_1 
is at most m l V , with equality if and only if the triangulation is unimodular (with respect to 
the lattice). 

In particular, the weighted efficiency of a triangulation of I 1 x A™ -1 is at most one, with 
equality for unimodular triangulations. 

Proof. If P is a lattice polytope (e.g., a cube), then the lattice spanned by tm is a sublattice 
of the one spanned by the point configuration P+ ■ ■ ■ +P, and coincides with it if and 
only if r is unimodular. Then, for unimodular triangulations the weighted size is then just 
the volume of P+ ■ ■ ■ +P, normalized to the unit parallelepiped, which equals m l V. For 
non- unimodular triangulations the weighted size is smaller than that. □ 

7. More on p ijm 
In this section we obtain the value of p2 7 m 

for any m and we will also show triangulations 
of J 3 x A 1 and I 3 x A 2 providing the values of ^3,2 and ^3^3 stated in Section [21 

7.1. Smallest weighted efficiency of triangulations of I 2 x A" 1-1 . Here we prove: 
Theorem 7.1. The smallest weighted efficiency p2.n1 of triangulations of I 2 x A" 1-1 is 

l\3m 2 /4] 

P2,m = \ 2 ' 
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That is, y/S/A for even m and -\/3/4 + 0(m 2 ) for odd m. A fine mixed subdivision of 
I 2 + ■ ■ ■ +I 2 corresponding to a triangulation of I 2 x A" 1-1 with that weighted efficiency is 
given in Figure^ 

Let Bm = B\+ ■ ■■ +B m be a cell in a fine mixed subdivision of I 2 + ™- +I 2 (a square of 
size m). Each Bi must be a simplex coming from the i-th copy of I 2 , and in order to have a 
fine mixed subdivision, the different B^s must lie in complementary affine subspaces. Then, 
there are the following possibilities: 

• Bm is a triangle, obtained as the sum of a triangle in one of the r's and a single 
point in the others. The weight of Bm is 1/2. 

• Bm is a quadrangle, obtained as the sum of two (non-parallel) segments from two 
of the BiS and a point in the rest of them. Three types of quadrangles can appear, 
depending on whether none, one or both of the two segments involved is a diagonal 
of J 2 : a square parallel to the axes, a rhomboid (both with area 1) and a diagonal 
square (of area 2). The weight of Bm is 1. 

In particular, the weighted size of the mixed subdivision equals T/2 + Si + S2 + R, 
where T, Si, S2 and R denote, respectively, the numbers of triangles, squares of area 1, 
squares of area 2 and rhombi in the subdivision. Since the total area of I 2 + ■ ■■ +I 2 is 
m 2 = T/2 + Si + 2S 2 + R, we conclude that: 

Proposition 7.2. The weighted size of a mixed subdivision of I 2 + ■ ■ ■ +I 2 equals m 2 — S2, 
where S2 is the number of squares of size 2 in the subdivision. 

With this we can already conclude that the fine mixed subdivisions shown in Figure 0] 
(one for m even and one for m odd) have weighted size equal to [3m 2 /4] , since they have 
exactly |_m 2 /4j squares of area 2. 



Figure 4. Fine mixed subdivisions corresponding to triangulations of J 2 x 
A m_1 with smallest weighted efficiency for m even (left) or odd (right). 

Our task is now to prove that no mixed subdivision can have more than |_m 2 /4j squares 
of area 2. For this we use: 

Lemma 7.3. Taking the i-th summands of all the mixed cells Bi+ ■ ■ ■ +B m in a mixed 
subdivision of Pi+ ■ ■ ■ +P m produces a polyhedral subdivision of Pi. If the mixed subdivision 
was fine, the polyhedral subdivision is a triangulation. 
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m 



Proof. By the Cayley Trick, the mixed subdivision of Pi+ ■ ■ ■ +P m induces a polyhedral 
subdivision of C(Pi, . . . , P m ), and a triangulation if the mixed subdivision is fine. Since the 
polytope Pi appears as a face in C(P\, . . . , P m ), any subdivision (resp. triangulation) of 
C(Pi, . . . , P m ) induces a subdivision (resp. triangulation) of Pi. That this subdivision is the 
one obtained taking the i-th summands of all the mixed cells follows from Theorem 16. II □ 

Proposition 7.4. Let Sm be a fine mixed subdivision of I 2 + ■ ■ ■ +I 2 . Let a and b be the 

number of summands L 2 's which are triangulated in one and the other possible triangulations 
of L 2 (i.e., using one or the other diagonal). Then, the number of squares of area 2 in Sm 
is at most ab and the weighted size of Sm is at least m 2 — ab. 

Proof. Each square of area 2 is the Minkowski sum of two opposite diagonals of two copies 
of I 2 (say, the ith and jth copies) , and a point in the other m — 2 of the copies of I 2 . Clearly, 
the two copies which contribute diagonals have to be triangulated in opposite ways. The 
only thing which remains to be shown is that the same pair of copies of I 2 cannot contribute 
two different squares of size 2. For this, observe that "contracting" in every mixed cell 
of a mixed subdivision all the summands other than the ith and jth should give a mixed 
subdivision of Pi + Pj . And in a mixed subdivision of two squares there is no room to put 
two different diagonal squares of area 2. □ 

In the statement of Proposition l7.4l we have that a + b = m. In particular, the maximum 
possible value of ab is [m 2 /A\. This finishes the proof of Theorem 17. II 

7.2. Smallest weighted efficiency triangulations of I 3 x A m , m = 2,3. Here we will 
try to visualize triangulations of L 3 x A 1 and I 3 x A 2 which minimize the weighted efficiency, 
which we computed using the integer programming approach sketched in Section |2 

Of course, we use the Cayley trick to decrease the dimension, so we show the corresponding 
fine mixed subdivision of a Minkowski sum instead of the triangulation itself. 



We have to give a subdivision of a 3-cube of size 2. For this we first cut the eight 
corners of the cube, producing a cubeoctahedron, a semi-regular 3-polytope with 6 
square and 8 triangular facets. The edges of the cubeoctahedron can be distributed 
in four "equatorial hexagons" each of which cuts the polytope into 2 halves. Figure 
depicts these two halves for one of the equatorial hexagons. The labels in the 
vertices are heights, interpreted as follows: Our point configuration is {0, 1, 2} 3 and 
the height of point k) is just i + j + k. 

It turns out that performing three of these four possible halvings, the cubeoctahe- 
dron is decomposed into six triangular prisms and two tetrahedra. In Figure El each 
half is actually decomposed into three prisms and one tetrahedron. These, together 
with the eight tetrahedra we have cut from corners, form a fine mixed subdivision 
of {0, 1, 2} 3 with 10 tetrahedra and 6 triangular prisms. Its weighted size is then 



• I 3 x A 1 



I 3 x A 1 = C(I 3 ,L 3 ) 



I 3 + 1 3 
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Figure 5. Fine mixed subdivision of I 3 + I 3 corresponding to a triangu- 
lation of I 3 x A 1 with smallest weighted efficiency. 



J 3 x A 2 | I 3 x A 2 = C(I 3 ,I 3 , 1 3 ) < — >I 3 + I 3 + I 3 

The fine mixed subdivision of I 3 + I 3 + I 3 is given in Figure El ft consists of 20 
triangular prisms, f6 tetrahedra and 2 parallelepipeds. 

Thus, the weighted size of the triangulation is 20^ + 16^ + 2y = 4p and then 
the smallest weighted efficiency is 



/44/3 




Figure 6. Fine mixed subdivision of I 3 + I 3 + I 3 corresponding to a tri- 
angulation of I 3 x A 2 with smallest weighted efficiency. 
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Let us explain how to interpret Figure Again, each point (i,j,k) in {0,1,2,3} 3 has 
been given a height i+j + k, which is written next to it. The subdivision is displayed in five 
parts. The left-top portion in the figure is a half cubeoctahedron exactly as the one in Figure 
El The reader has to assume it divided into a tetrahedron and three triangular prisms, as 
before. The left-bottom portion consists of a corner tetrahedron and three triangular prisms 
located at corners of the cube, each joined to the half cubeoctahedron by a tetrahedron. 
So far we have six prisms and five tetrahedra, and the same is got from the right-top and 
right-bottom portions of the Figure. 

In between the two half cubeoctahedra, however, we have now a hexagonal prism de- 
composed into two triangular prisms and two quadrilateral prisms, all of height The 
hexagonal prism is surrounded by a belt formed by six triangular prisms and six tetrahedra. 

The thick edges in Figure El represent edges of the big 3-cube of side 3, whose three visible 
facets are drawn by dotted lines. We have also shaded the facets of the subdivision which 
are contained in those facets of the big 3-cube. 
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